Labyrinthine pathways towards supercycle attractors in unimodal maps 
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We uncover previously unknown properties of the family of periodic superstable cycles in uni- 
modal maps characterized each by a Lyapunov exponent that diverges to minus infinity. Amongst 
the main novel properties are the following: i) The basins of attraction for the phases of the cycles 
develop fractal boundaries of increasing complexity as the period-doubling structure advances to- 
wards the transition to chaos, ii) The fractal boundaries, formed by the preimages of the repellor, 
display hierarchical structures organized according to exponential clusterings that manifest in the 
dynamics as sensitivity to the final state and transient chaos, iii) There is a functional composition 
renormalization group (RG) fixed-point map associated to the family of supercycles. iv) This map 
is given in closed form by the same kind of g-exponential function found for both the pitchfork and 
tangent bifurcation attractors. v) There is a final stage ultra-fast dynamics towards the attractor 
with a sensitivity to initial conditions that decreases as an exponential of an exponential of time. 
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I. INTRODUCTION 

The family of superstable orbits - often called super- 
cycles - of unimodal maps have played a prominent role 
in the historical uncovering of the cascade of period- 
doubling bifurcations that leads to the transition to chaos 
with the celebrated universal features displayed by the 
Feigenbaum attractor 0, [I]- These periodic attrac- 
tors are easy to find because one of its cycle positions 
corresponds to the maximum of the map. And then, 
they are convenient for basic computational purposes be- 
cause trajectories originating at different positions settle 
very fast into the attractor for the reason that their Lya- 
punov exponent diverges to minus infinity. And so they 
offered a suitable path for the initial characterization of 
the period-doubling route to chaos. Oddly enough, sev- 
eral features associated to the dynamical evolution to- 
ward these attractors have remained unexplored. Here 
we provide an account of these properties, some shared 
by other families of periodic attractors and some specific 
of the supercycles. We discuss explicitly the case of a 
map with quadratic maximum but the results are easily 
extended to general nonlinear! ty z > 1. In closing we 
indicate the broader theoretical interest in these proper- 
ties. 

In order to obtain dynamical properties with previ- 
ously unstated detail we determined the organization of 
the entire set of trajectories as generated by all possi- 
ble initial conditions. We find that the paths taken by 
the full set of trajectories in their way to the supercy- 
cle attractors (or to their complementary repellors) are 
far from unstructured. The preimages of the attractor of 
period 2^, N = 1,2,3,... are distributed into different 
basins of attraction, one for each of the 2^ phases (posi- 
tions) that compose the cycle. When N > 2 these basins 



are separated by fractal boundaries whose complexity in- 
creases with increasing TV. The boundaries consist of the 
preimages of the corresponding repellor and their posi- 
tions cluster around the 2^ — 1 repellor positions accord- 
ing to an exponential law. As N increases the structure 
of the basin boundaries becomes more involved. Namely, 
the boundaries for the 2^ cycle develops new features 
around those of the previous 2^^^ cycle boundaries, with 
the outcome that a hierarchical structure arises, leading 
to embedded clusters of clusters of boundary positions, 
and so forth. 

The dynamics associated to families of trajectories al- 
ways displays a distinctively concerted order that reflects 
the repellor preimage boundary structure of the basins 
of attraction. That is, each trajectory has an initial 
condition that is identified as an attractor (or repellor) 
preimage of a given order, and this trajectory necessarily 
follows the steps of other trajectories with initial condi- 
tions of lower preimage order belonging to a given chain 
or pathway to the attractor (or repellor). This feature 
gives rise to transient chaotic behavior different from that 
observed at the last stage of approach to the attractor. 
When the period 2^ of the cycle increases the dynamics 
becomes more involved with increasingly more complex 
stages that reflect the preimage hierarchical structure. 
As a final point, at the closing part of the last leg of the 
trajectories an ultra-rapid convergence to the attractor 
is observed, with a sensitivity to initial conditions that 
decreases as an exponential of an exponential of time. In 
relation to this we find that there is a functional com- 
position renormalization group (RG) fixed-point map as- 
sociated to the supercycle attractor, and this can be ex- 
pressed in closed form by the same kind of q -exponential 
function found for both the pitchfork and tangent bifur- 
cation attractors @, like that originally derived by 



2 



Hu and Rudnick for the latter case [y] . 

Before proceeding to give details in the following sec- 
tions of the aforementioned dynamics we recall basic 
features of the bifurcation forks that form the period- 
doubling cascade sequence in unimodal maps, epitomized 
by the logistic map Li{x) = 1 — ^x^, —l<x<l, 
O</i<2[l|,0,[3]. The superstable periodic or- 
bits of lengths 2^, N = 1,2,3,..., are located along 
the bifurcation forks, i.e. the control parameter value 
M — 'Pn < Moo for the superstable 2^-attractor is that 
for which the orbit of period 2^ contains the point a; = 0, 
where fioo — 1.401155189... is the value of fi at the 
period-doubling accumulation point. The positions (or 

phases) of the 2^-attractor are given by x„i — /1™'(0), 

m — 1,2,..., 2^. Notice that infinitely many other se- 
quences of superstable attractors appear at the period- 
doubling cascades within the windows of periodic attrac- 
tors for values of /i > . Associated to the 2^-attractor 
at /X = Jij^ there is a (2^ — l)-repellor consisting of 
2^ — 1 positions j/„j, m = 1, 2, . . . , 2^ — 1. These po- 



sitions are the unstable solutions, 
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, of 2/ = >iy), n = 1,2,. ..,7V. The first, n - 1, 

originates at the initial period-doubling bifurcation, the 
next two, n = 2, start at the second bifurcation, and so 
on, with the last group of 2^^^, n — N, stemming from 
the — 1 bifurcation. See Fig. [TJ We find it useful to 
order the repellor positions, or simply, repellors, present 
at /i = /Iat, according to a hierarchy or tree, the 'oldest' 
with n = 1 up to the most recent ones with n = N . The 
repellors' order is given by the value of n. Finally, we de- 
fine the preimage x^^^ of order k of position x to satisfy 
X — /i*^'^)(x^*'^) where h''^\x) is the fc-th composition of 

the map h(x) = /A.^ {^)- We have omitted reference 
to the 2^-cycle in x'^^'^ to simplify the notation. 



II. PREIMAGE STRUCTURE OF 
SUPERCYCLE ATTRACTORS 

The core source of our description of the dynamics to- 
wards the supercycle attractors is a measure of the rel- 
ative 'time of flight' tf{xo) for a trajectory with initial 
condition xg to reach the attractor. The function tf{xo) 
is obtained for an ensemble representative of all initial 
conditions — 1 < xq < 1. This comprehensive informa- 
tion is determined through the numerical realization of 
every trajectory, up to a small cutoff £ > at its final 
stage. The cutoff e regards a position Xf < Xm ± e to 
be effectively the attractor phase Xm- This, of course, 
introduces an approximation to the real time of flight, 
that can be arbitrarily large for those xq close to a re- 
pellor position Urn or close to any of its infinitely many 

(k) 

preimages, xin , k — 1,2,..., N > 2. In such cases 
the finite time tf{xo;e) can be seen to diverge — > oo 
as xq Um and e — s- 0. As a simple illustration, in 
Fig. [5] we show the time of fiight tf{xQ) for the period 
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FIG. 1: Solutions of a; = /Jf \x) , N ^ 0, 1, 2, where 
is the logistic map, as a function of the control parameter 
/i. The attractors (stable solutions) become unstable at each 
bifurcation, turning into repellors. When for a given value 
of fi the attractor has 2^ positions there are 2^^ — 1 repellor 
positions. 



2 supercycle at /z = /x^ with e = 10 ^ together with 
the twice-composed map f-^^{x) and a few representa- 
tive trajectories. We observe two peaks in tf{xQ]e) at 
2/1 = (-1 -I- VI + 47Ii)/27ri ~ 0.6180340..., the fixed- 
point repellor and at its (only) preimage xj^'' 



where yi = /^^''(x^^''). Clearly, there are two basins of 
attraction each for the two positions or phases, xi — 
and X2 = 1, of the attractor. For the former it is the 
interval x']^^ < xq < yi, whereas for the second one con- 
sists of the intervals — 1 < xo < x^^'' and yi < xq < 1. 
The multiple-step structure of tf{xo), of four time units 
each step, reflects the occurrence of intervals of initial 
conditions with common attractor phase preimage order 
k. 

For the next supercycle - period 4 - the preimage struc- 
ture turns out to be a good deal more involved than the 
straightforward structure for /Ii • In Fig. [3] we show the 
times of flight tf{xo) for the N — 2 supercycle at /j, = 
with e = 10"^, the map f^^{x) is superposed as a ref- 
erence to indicate the four phases of the attractor (at 
xi = 0, X2 = 1, X3 ~ -0.3107 . . ., and x^ ~ 0.8734 . . .) 
and the three repellor positions (at yi ~ 0.5716635 . . ., 
2/2 ~ 0.952771..., and 2/3 - -0.189822...). In Fig. [3] 
there are also shown four trajectories each of which ter- 
minates at a different attractor phase. We observe a pro- 
liferation of peaks and valleys in ^^(xo), actually, an infl- 
nite number of them, that cluster around the repellor at 
2/1 ~ 0.5716635 . . . and also at its preimage at xj^^ = — j/i 
(these are the positions at fi = 'p2 of the 'old' repellor 
and its preimage in the previous A'^ = 1 case). Notice 
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FIG. 2: Left axis: The twice iterated map fji {x), 1 
(in grey line). Right axis: Time of flight tf{x), the number 
of iterations necessary for a trajectory with initial condition 
at X to reach an attractor position. The values of x near the 
high spikes correspond to initial conditions very close to the 
repellor and its preimage. We present three example trajec- 
tories (and the y = x line as an aid to visualize them): The 
dotted line shows a trajectory that starts at the attractor po- 
sition X = and remains there. The solid line is a trajectory 
starting near the repellor at yi, and after a large number of 
iterations reaches the attractor position x — 1. Finally, the 
dash-dotted line is an orbit starting at a; = 0.5 that in just a 
few iterations reaches x — 0. 



that the steps in the valleys of tf{xQ) are now eight time 
units each. The nature of the clustering of peaks (re- 
pellor phase preimages) and the bases of the valleys (at- 
tractor phase preimages) is revealed in Fig. |4] where we 
plot tf in a, logarithmic scale for the variables ±(x — yi). 
There is an exponential clustering of the preimage struc- 
ture around both the old repellor and its preimage. This 
scaling property is corroborated in Fig. [5] from where we 
obtain x-yi ~ 7.5x 10"^ exp(0.80 0, where ^ = 1,2,3, .. . 
is a label for consecutive repellor preimages. 

A comparable leap in the complexity of the preimage 
structure is observed for the following - period 8 - super- 
cycle. In Fig. [6] we show tf{xo) for the iV = 3 super- 
cycle at /J, = /I3 with £ — 10^^, together with the map 
fj^^ (x) placed as reference to facilitate the identification 
of the locations of the eight phases of the attractor, xi 
to Xg, and the seven repellor positions, yi to j/y. Besides 
a huge proliferation of peaks and valleys in tf{xo), we 
observe now the development of clusters of clusters of 
peaks centered around the repellor at yi ~ 0.5626447 . . . 
and its preimage at a;^^'' = —yi (these are now the po- 
sitions of the original repellor and its preimage for the 
N — 1 case when ^ = 7I3). The steps in the valleys 
of tf(xQ) have become sixteen time units each. Simi- 
larly to the clustering of peaks (repellor phase preim- 
ages) and valleys (attractor phase preimages) for the pre- 




FIG. 3: Same as Fig. [2] but for fjij{x), JL2 ^ 1.31070264. 
Two orbits (in solid bold and dashed bold lines) start very 
near each other and by the position j/i of the old repellor, 
indistinguishable at the beginning they take very different 
subsequent paths to reach the attractor positions at a; = 1 
and X ~ 0.8734. See bottom right panel. Another two orbits 
(in dotted and solid lines) start at an attractor position x = 1 
and at a repellor preimage positions, respectively. Finally, one 
more orbit (in dashed line) starts at an intermediate initial 
condition an reaches very quickly the attractor position at 
X ~ —0.310703. See bottom left panel. 



vious supercycle at 7I2, the spacing arrangement of the 
new clusters of clusters of peaks is determined in Fig. [7| 
where we plot t/ in a logarithmic scale for the variables 
±(x— t/i). In parallel to the previous cycle an exponential 
clustering of clusters of the preimage structure is found 
around both the old repellor and its preimage. This scal- 
ing property is quantified in Fig. [8] from where we obtain 
x-yi ~ 8.8xl0"^exp(0.84?), where I = 1,2,3, . . . counts 
consecutive clusters. 

An investigation of the preimage structure for the next 
= 4 supercycle at ^ — JI^ leads to another substantial 
increment in the involvedness of the structure of preim- 
ages but with such density that makes it cumbersome 
to describe here. Nevertheless it is clear that the main 
characteristic in the dynamics is the development of a hi- 
erarchical organization of the preimage structure as the 
period 2^ of the supercycles increases. 



III. FINAL STATE SENSITIVITY AND 
TRANSIENT CHAOS 

With the knowledge gained about the features dis- 
played by the times of flight tf{xo) for the first few super- 
cycles it is possible to determine and understand how the 
leading properties in the dynamics of approach to these 
attractors arise. The information contained in tf{xo) can 
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FIG. 4: Top panel: Time of flight tf{x) for TV = 2 (as in Fig. 
3), the black lines correspond to initial conditions that ter- 
minate at the attractor positions a; = and x ~ —0.310703, 
while the grey lines to trajectories ending at a:: = 1 and 
X ~ 0.8734. Right (left) bottom panel: Same as top panel, 
but plotted against the logarithm of x — (yi—x). It is ev- 
ident that the peaks are arranged exponentially around the 
old repellor position yi, i.e., they appear equidistant in a log- 
arithmic scale. See Fig. [S] 
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FIG. 5: Corroboration of preimage exponential clustering 
around the repellor position yi ~ 0.571663 when N = 2 and 
7I2 — 1.31070264. The variable I labels consecutive equidistant 
peaks in the inset. The peaks correspond to preimages of the 
repellor. See text. 



be used to demonstrate in detail how the concepts of fi- 
nal state sensitivity 01 - due to attractor multiplicity 
- and transient chaos [4I - prevalent in the presence of 
repellors that coexist with periodic attractors - realize 
in a given dynamics. Final state sensitivity is the con- 




FIG. 6: Same as Fig. |3]for fj, 'ix), Jl^~ 1.38154748. Here 
yi= 0.56264475. ^ 




FIG. 7: Same as Fig. H for = 3. The black lines corre- 
spond to initial conditions that terminate at any of the four 
attractor positions close or equal to a; = 0, while the grey lines 
to trajectories ending at any of the other four attractor po- 
sitions close or equal to a; = 1. As the bottom panels show, 
in logarithmic scale, in this case there are (infinitely) many 
clusters of peaks (repellor preimages) equidistant from each 
other. 



sequence of fractal boundaries separating coexisting at- 
tractors. In our case there is always a single attractor 
but its positions or phases play an equivalent role 
Transient chaos [2] is due to fast separation in time of 
nearby trajectories by the action of a repellor and results 
in a sensitivity to initial conditions that grows exponen- 
tially up to a crossover time after which the decay sets 
in. We describe below how both properties result from an 
extremely ordered flow of trajectories towards the attrac- 
tor. This order is imprinted by the preimage structure 



5 




FIG. 8: Same as Fig. E] for iV = 3. 
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FIG. 9: Top panel: Detail of the map /l^ {x) and a tra- 
jectory in the proximity of the attractor position located at 
X ~ 0.850780. The trajectory originated very close to the re- 
pellor position at a; = 0.819378. Bottom panel: The same 
trajectory of the eight-iterated map as a function of time t. 
Inset: Corroboration of the exponential nature of the trajec- 
tory after it leaves the repellor and before the final approach 
to the attractor position. 



described in the previous section. 

For the simplest supercycle at /i = there is trivial 
final state sensitivity as the boundary between the two 
basins of the phases, xi = and 2:2 = 1, consists only of 
the two positions, yi ~ 0.5716635 . . . and x^i^ — —yi See 
Fig. [2I Consider the length 5 of a small interval around a 
given value of xq containing either yi or xj^', when (5 — s- 
any uncertainty as to the final phase of the trajectory 
disappears. It is also simple to verify that when xq is 



close to yi or the resulting trajectories increase their 
separation at initial and intermediate times displaying 
transient chaos in a straightforward fashion. In Fig. [2]we 
show the same type of transitory exponential sensitivity 
to initial conditions for a trajectory at ^ = JI^ after it 
reaches a repellor position in the final journey towards 
the period 8 attractor. This behavior is common to all 
periodic attractors when trajectories come near a repellor 
at the final leg of their journey. 




FIG. 10: Time of flight t j {x) ior N = 2 (equivalent to Fig. Ell . 
Different types of lines are used for different sub-basins of at- 
traction that form the fractal boundary between the attractor 
positions. The black solid line corresponds to the attractor 
position a; = 0, the grey solid line to a; = 1, dashed grey line to 
X ~ -0.310703, and dashed black line to a; ~ 0.873470. The 
sub-basins are separated in an alternated fashion by repulsor 
preimages (peaks) and there are just two types of sub-basins 
depending on |a:| being larger or smaller than y\. 

For II = 'P2 there are more remarkable properties aris- 
ing from the more complex preimage structure. There 
is a concerted migration of initial conditions seeping 
through the boundaries between the four basins of at- 
traction of the phases. These boundaries, shown in Fig. 
[TOl form a fractal network of interlaced sub-basins sepa- 
rated from each other by two preimages of different re- 
pellor phases and have at their bottom a preimage of an 
attractor phase. Trajectories on one of these sub-basins 
move to the nearest sub-basin of its type (next-nearest 
neighbor in actual distance in Fig. [10]) at each iteration of 

the map fj^^ (4 time steps for the original map fji^ ) . The 
movement is always away from the center of the cluster 
at the old repellor position j/i or at its preimage x^^^* (lo- 
cated at the steepest slope inflection points of f-jj^ shown 
in Fig. Once a trajectory is out of the cluster (con- 
tained between the maxima and minima of f— next to 
the mentioned inflection points) it proceeds to the basin 
of attraction of an attractor phase (separated from the 
cluster by the inflection points with gentler slope of f-^^ 
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in Fig. [3]) where its final stage takes place. When con- 
sidering a large ensemble of initial positions, distributed 
uniformly along all phase space, the common journey to- 
wards the attractor displays an exceedingly ordered pat- 
tern. Each initial position xq within either of the two 
clusters of sub-basins is a preimage of a given order k of 
a position in the main basin of attraction. Each iteration 
of fj^^ reduces the order of the preimage from k to k — 1, 

and the new position Xq — f^^(xo) replaces the initial 
position Xq (a preimage of order fc — 1) of another trajec- 
tory that under the same time step has migrated to the 
initial position Xq (a preimage of order A; — 2) of another 
trajectory, and so on. 



120- 



100- 



/eo- 




FIG. 11; Evolution of an orbit starting very close to a 
peak (repellor preimage) inside a cluster when A*' = 2, /Ij ~ 
1.31070264 and yi= 0.571663. (See Fig. Ell .The trajectory 
describes an exponential in time moving away from yi. In 
the figure, positions of the trajectory correspond to the black 
solid dots, and iterations correspond to the number associated 
to each dot. 

It is clear that the dynamics at ii — 'p2 displays sen- 
sitivity to the final state when the initial condition xq is 
located near the core of any of the two clusters of sub- 
basins at yi and x^^^ that form the boundary of the at- 
tractor phases. Any uncertainty on the location of xq 
when arbitrarily close to these positions implies uncer- 
tainty about the final phase of a trajectory. There is also 
transient chaotic behavior associated to the migration of 
trajectories out of the cluster as a result of the organized 
preimage resettlement mentioned above. Indeed, the ex- 
ponential disposition of repellor preimages shown in Fig. 
[5] is actually a realization of two trajectories with ini- 
tial conditions in consecutive peaks of the cluster struc- 
ture. Therefore, the exponential expression given in the 
previous Section in relation to Fig. [5] can be rewritten 
as the expression for a trajectory, Xr — exp(Ae//T), 
with, Xr — X — yi, Xq ~ 7.5 X 10~^, and Ae// = 6.4, 
where r = 1,2,3, .. .. Straightforward differentiation of 
Xr with respect to xq yields an exponential sensitivity to 



initial conditions with positive effective Lyapunov coeffi- 
cient Ae//. 

As it can be anticipated, the dynamics of approach 
to the next supercycle at = 7I3 can be explained by 
enlarging the description presented above for ^ = /I2 
with the additional features of its preimage structure al- 
ready detailed in the previous section. As in the previous 
case, trajectories with initial conditions xq located inside 
a cluster of sub-basins of the attractor phases will pro- 
ceed to move out of it in the systematic manner described 
for the only two isolated clusters present when /i = /Ij. 
Only now there is an infinite number of such clusters ar- 
ranged into two bunches that group exponentially around 
the old repellor position yi and around its preimage x^i'^ . 
See Figs. [6] to HI Once such trajectories leave the cluster 
under consideration they enter into a neighboring cluster, 
and so forth, so that the trajectories advance out of these 
fractal boundaries through the prolonged process of mi- 
gration out of the cluster of clusters before they proceed 
to the basins of the attraction of the eight phases of this 
cycle. In Fig. [T^we show one such trajectory in consec- 
utive times t — 1, 2, 3, . . . for the original map and also 
in multiples of time t — 2^, 22^, 32^, . . . The logarithmic 
scale of the figure makes evident the retardation of each 
stage in the process. As when ^ = JI2, it is clear that in 
the approach to the 11 ='p^ attractor there is sensitivity 
to the final state and transitory chaotic sensitivity to ini- 
tial conditions. Again, the exponential expression, given 
in the previous Section associated to the preimage struc- 
ture of clusters of clusters of sub-basins shown in Fig. 
[21 can be interpreted as the expression of a trajectory 
of the form Xr — exp(Ae//T). Differentiation of Xr 
with respect to xq yields again an exponential sensitiv- 
ity to initial conditions with positive effective Lyapunov 
coefficient Ae//. 

When the period 2^ of the cycles increases we observe 
that the main characteristic in the dynamics is the de- 
velopment of a hierarchical organization in the flow of an 
ensemble of trajectories out of an increasingly more com- 
plex disposition of the preimages of the attractor phases. 
The role played by a cluster of sub-basins when n = JI2 
is preserved for the clusters present when fj, — Jl^ but the 
additional presence of clusters of clusters of sub-basins 
introduces a similar role in the dynamics at a broader 
scale in the preimage structure. A complication that will 
arise in a comparable fashion every time N increases. 



IV. SUPER STRONG INSENSITIVITY TO 
INITIAL CONDITIONS 

We turn now to discuss the closing point of the ap- 
proach to the supercycle attractors and find out what is 
the form of the sensitivity to initial conditions in the 
t 3> 1 limit, as required in its usual definition. The ques- 
tion is pertinent because for these attractors the ordinary 
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FIG. 12: A trajectory for Jl^~ 1.38154748. Top panel: The 
circles are positions for consecutive times in the iterations of 
the map fji.^{x). The orbit starts very close to the old (period 
one) repellor j/i , then moves close to a period two and subse- 
cuently to a period four repellor, until finally the trajectory 
arrives at the period eight attractor. Bottom panel: Selec- 
tion of a time subsequence of multiples of 2^ shows the same 
trajectory as an evolution towards one particular attractor 
position. 



diverges to minus infinity {dfj^j^{x)/dx = at a; = 0) and 
cannot have an exponential form — exp(Af) with 
A < 0. 

Representative resuhs for the last segment of a trajec- 
tory and the corresponding sensitivity obtained from 
a numerical investigation are shovsfn in Fig. 1131 Only the 
first two steps of a trajectory with xq = 0.1 of the map 
/^^ ^ can be seen in Fig. [13] (a). A considerable enlarge- 
ment of the spatial scale (that requires computations of 
extreme precision [9]) makes it possible to observe a to- 
tal of seven steps (56 iterations in the original map), as 
shown with the help of logarithmic scales in Fig. [13] (b). 
Fig. [13] (c) shows both the same trajectory and the sen- 
sitivity S,t = dxt/dxQ in a logarithmic scale for xt and 
£,t and a normal scale for the time t. The trajectory is 
accurately reproduced (indistinguishable from the curve 
in Fig. [13] (c)) by the expression xt = u^^ exp(6expci), 
b = Inuxo (with xq > 0) and c ~ In 2, where u > 

is obtained from the form I fl^ I ~ ux'^ that the 2^-th 
composed map takes close to x — 0. This expression 
for Xt is just another form of writing uxt — (uxq)^ , the 
result of repeated iteration of mx^. For the logarithm 
of the sensitivity we have In — — Ihxq + t\Ti2 + hixt 
where the last (large negative) term dominates the first 
two. Thus, we find that the sensitivity decreases faster 
than an exponential, and more precisely, decreases as the 
exponential of an exponential. 

It is interesting to note that, associated to the gen- 
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FIG. 13: Panel (a): Detail of a trajectory in its final stage of 
approach to the a; = attractor position (dotted line) for the 

map /A. ' (x) (solid line). In this example A'' = 3. Panel (b): 
Same as (a) in double logarithmic scale. Panel (c): Trajectory 
Xt and sensitivity to initial conditions in logarithmic scale 
versus time t. Both functions are indistinguishable. Note that 
there is an ultra-rapid convergence, of only a few time steps, 
to the origin a: = 0. See text. 



borhood of a; = 0, there is a map f*{x) that satis- 
fies the functional composition and rescaling equation 
f*{f*{x)) = a~'^ f*{ax) for some finite value of a and 
such that f*{x) — ux"^ + o{x^). The fixed-point map 
f*{x) possesses properties common to all superstable at- 
tractors of unimodal maps with a quadratic extremum. 
Indeed, there is a closed form expression that satisfies 
these conditions. This is 

f*{x) = xexpg{u''~'^x), 

where exp^(a;) is the ^-exponential function exp^(a;) = 
[1 -I- (1 — q)x]^^^^~''^K The fixed-point map equation is 
satisfied with a — 2^/^'^^) and q = 1/2. The same type 
of RG solution has been previously found to exist for 
the tangent and pitchfork bifurcations of unimodal maps 
with general nonlinearity > 1 0], [l], Use of the 
map f*{x) reproduces the trajectory Xt/xo and sensibil- 
ity ^t shown in Fig. [13] (c). 



V. SUMMARY 

We studied the properties of the first few members of 
the family of superstable attractors of unimodal maps 
with quadratic maximum and obtained a precise under- 
standing of the complex labyrinthine dynamics that de- 
velops as their period 2^ increases. The study is based 
on the determination of the function tf{xo), the time of 
flight for a trajectory with initial condition xq to reach 
the attractor or repellor. The novelty and scope of the 
study is that the function t/(a^o) was determined for all 
initial conditions a;o in a partition of the total phase space 
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— 1 < a;o ^ Ij and this provides a complete picture for 
each attractor-repellor pair. We observed how the frac- 
tal features of the boundaries between the basins of at- 
traction of the positions of the periodic orbits develop 
a structure with hierarchy, and how this in turn reflects 
on the properties of the trajectories. The set of trajec- 
tories produce an ordered flow towards the attractor or 
towards the repellor that reflect the ladder structure of 
the sub-basins that constitute the mentioned boundaries. 
As 2^^ increases there is sensitivity to the final position 
for almost all Xq, and there is a transient exponentially- 
increasing sensitivity to initial conditions for almost all 
xq. We observed that transient chaos is the manifesta- 
tion of the trajectories' controlled flow out of the fractal 
boundaries, that suggest that for large 2^ the flow be- 
comes an approximately self-similar sequence of stages. 
Finally, we looked at the closing segment of trajectories 
at which a very fast convergence to the attractor posi- 
tions occurs. Here we found 'universality class' features, 
as the trajectories and sensitivity to initial conditions 
are replicated by an RG fixed-point map obtained un- 
der functional composition and rescaling. This map has 
the same g-deformed exponential closed form found to 



hold also for the pitchfork and tangent bifurcations of 
unimodal maps 

The new knowledge gained offers a significant qualita- 
tive and quantitative increment in the understanding of 
the dynamics at the transition to chaos, in this case via 
the period doubling route. The dynamics of approach to 
the Feigenbaum attractor and repellor, the supercycle of 
period 2°°, is, to our knowledge, not known at the level 
of detail given here for the first few supercycles. These 
comprehensive properties of the Feigenbaum attractor 
and repellor are presented [isj in tandem with this 
work. There is recent renewed interest in the dynamics 
of low-dimensional nonlinear maps as access to exact 
and detailed information for specific model systems is 
important in assessing theoretical developments, such as 
generalizations of ordinary statistics. There are positive 
indications that the multifractal critical attractors 
present in these maps play such role [l^], [IH, [13] ■ 

Acknowledgments. D.S. and A.R. are grateful to 

Hugo Hernandez Saldana for interesting discussions. 

Partial support by DGAPA-UNAM and CONACyT 
(Mexican Agencies) is acknowledged. 



[1] H.G. Schuster, Deterministic Chaos. An Introduction, 
2nd Revised Edition, VCH Publishers, Weinheim, 1988. 

[2] C. Beck, F. Schlogl, Thermodynamics of Chaotic Sys- 
tems, Cambridge University Press, UK, 1993. 

[3] R.C. Hilborn, Chaos and Nonlinear Dynamics, 2nd Re- 
vised Edition, Oxford University Press, New York, 2000. 

[4] A. Robledo, Physica A 314, 437 (2002); A. Robledo, 
Physica D 193, 153 (2004). 

[5] F. Baldovin, A. Robledo, Europhys. Lett. 60, 518 (2002). 

[6] B. Hu and J. Rudnick, Phys. Rev. Lett. 48, 1645 (1982). 

[7] C. Grebogi, S.W. McDonald, E. Ott, J.A. Yorke, Phys. 
Lett. 99A, 415 (1983). 

[8] S.W. McDonald, C. Grebogi, E. Ott, J.A. Yorke, Physica 



17D, 125 (1985). 
[9] Our simulations used MAPM, the 

arbitrary precision C-|— |- library 

(|http://www.tc.umn.edu^ringx004/mapm-main.html ) . 

We implemented our calculations taking between 200 

and 1200 significant digits. 
[10] E. Mayoral, A. Robledo, Phys. Rev. E 72, 026209 (2005). 
[11] H. Hernandez-Saldana, A. Robledo, Physica A 370, 286 

(2006). 

[12] A. Robledo, Physica A 370, 449 (2006). 
[13] A. Robledo, L.G. Moyano, Phys. Rev. E (submitted); 
/cond-mat/XXXX. 



